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C^ ! Abstract 

The parity type of a Latin square is defined in terms of the numbers of even and 

, odd rows and columns. It is related to an Alon- Tarsi- like conjecture that applies 

^ I to Latin squares of odd order. Parity types are used to derive upper bounds for 

cn ' the size of autotopy groups. A new algorithm for finding the autotopy group of a 

^ . Latin square, based on the cycle decomposition of its rows, is presented, and upper 

Q I bounds for the size of autotopy groups are derived from it. 

Q . Keywords: Latin square, parity type, cycle structure, autotopy group. 
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1 Introduction 



H ' For a given positive integer n let [n] denote the set {1, . . . , n}. A Latin square of order n 

is an n X n array of numbers in [n] so that each number appears exactly once in each row 
and column. Let LS{n) be the set of Latin squares of order n and let Sn be the symmetric 
group of permutations of [n]. The group 5^ = ^n x S^ x Sn acts on LS{n) by isotopism. 
An isotopism is a triple (a,/3,7) G S^ that acts on LS{n), so that a permutes the rows, 
P permutes the columns and 7 permutes the symbols of a given Latin square. Two Latin 
squares are called isotopic if there is an isotopism that transforms one to the other. Being 
isotopic is an equivalence relation and the set of Latin squares that are isotopic to each 
other is an isotopy class. Let J^„ = 5"^ be the group of isotopisms on LS{n), ^„,(L) 
the isotopy class of L, and 2l(L) = {0 G J^„|6(L) = L} the autotopy group of L (also 



*Thc author thanks Ian Wanless for providing a large part of the code needed to perform the computer 
simulations, and an anonymous referee for many helpful comments. 
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denoted in the literature by Is(L) [12], Atp(L) [2, 15], and Atop(L) [13]), whose members 
are called autotopisms of L. We have that |c/„| = (?i!)'^ and |^„(L)| = (n!)'^/ |2l(L)|. 
The rows and columns of a Latin square can be viewed as permutations of the ordered 
set (1, 2, . . . , n). A reduced (also called normalized) Latin square is a Latin square whose 
first row and first column are the identity permutation. Let RLS{n) be the set of reduced 
Latin squares of order n. For L e RLS{n) let J^^^ = {0 G J^„|6(L) G RLS{n)} and 
let ^ni-^) — '^ni.L) n RLS{n). Any Latin square can be transformed by isotopism to a 
reduced Latin square by permuting the rows to form the identity permutation in the first 
column and then permuting the columns 2, . . . , n to form the identity permutation in the 
first row. Thus, each reduced Latin square is isotopic to n\{n — 1)! distinct Latin squares 
in this way. It follows that |J^„(L)| / \^^{L)\ = n\(n — 1)! and hence \^ni\ = n ■ n\ and 

|j^^(L)|=n-n!/|2l(L)| (1.1) 

For these and other facts about isotopisms and autotopisms the reader is referred to Denes 
and A. D. Keedwell's book [4] and to the works of Janssen [10], McKay and Wanless [12] 
and McKay, Meynert and Myrvold [11] among many others. 

The sign (or parity) function sgn : Sn — ?■ { — 1,1} is defined so that sgn(o") = 1 if and 
only if a is an even permutation. The parity of a Latin square L, denoted par(L), is the 
product of the signs of all its rows and columns. A Latin square is even (resp. odd) if its 
parity is 1 (resp. -1). Let ELS{n) (resp. OLS{n)) be the set of even (resp. odd) Latin 
squares of order n. Let RELS{n) and ROLS{n) be the corresponding sets of reduced 
Latin squares. If n is odd \ELS{n)\ = \OLS{n)\. If n is even, a known conjecture of 
Alon and Tarsi [1] states that \ELS{n)\ ^ \OLS{n)\. It was also conjectured in [15] 
that \RELS{n)\ ^ \ROLS{n)\ for all n (for even n this is equivalent to the Alon- Tarsi 
conjecture). For further results regarding even and odd Latin squares, the reader is 
referred to the works of Janssen [10], Drisko [5, 6], Zappa [17], Glynn [8] and Stones and 
Wanless [15]. 

This work is mainly concerned with the relation between parities of Latin squares and of 
their rows and columns and autotopy groups. In Section 2 an expression for the members 
of J^^ ^ is derived and the notions of parity type and parity class are introduced. These 
lead to expressions for the numbers of even and odd reduced Latin squares of odd order 
in an isotopy class. These expressions, although stated and proved differently coincide 
with results of Stones and Wanless [15]. In Section 3 an extension of the Alon- Tarsi 
Latin square conjecture is discussed and a related conjecture on parity types is presented. 
In Sections 4 and 5 bounds are obtained for the order of the autotopy group 2t(L). In 
Section 4 the bounds are obtained using parities of the rows and columns of L and in 
Section 5 the bounds are derived from the cycle structure of the rows of L. 
Convention: Throughout the manuscript, when viewing rows and columns of a Latin 
square as permutations in Sn, the following convention will hold: the number i appearing 
in the jth place of the row (column) a of L signifies that a{i) = j. By this convention, in 
order to transform a column (or row) a to the identity permutation, we have to apply the 
permutation a~^ to the rows (resp. columns) of L (see Lemma l(i)). It is also assumed 
that permutations are applied from right to left. 
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2 Parity types and parity classes 

A permutation acting on a Latin square acts on another permutation (a row or a column) 
in two different ways, as described by the following lemma. Recall the convention of this 
paper that when i appears in the jth place of a row (or column) ir it means that 7r{i) = j. 

Lemma 1. Let L be a Latin square of order n. 

(i) The result of permuting the rows (resp. columns) of L, hy a E Sn, on any column 
(resp. row) vr is an. 

(ii) The result of permuting the symbols of L, by a E Sn, on any row or column vr is 
7ra~^ . 

Proof, (i) Let tt be a row or a column. Suppose 7r{i) = j and a{j) = k then i appears in 
the jth place of vr, and after applying a, the number in the jth place of vr (i in this case) 
moves to the kth place. Thus the resulting permutation takes i to k, so it is air. 
(ii) Denote by a(vr) the permutation obtained by applying a to the each symbol In the 
row (or column) vr. Suppose 7i{i) = j and a{i) = k. Since i is in the jth place of vr, k is 
in the jth place of a(vr). Thus a(vr)(/c) = j. It follows that a(vr) = 7ia~^. D 

The following proposition appears in a different form as part of the proof of Theorem 
2.1 in [15]. It describes the n ■ n\ elements in the set ^^ l- 

Proposition 1. Let L be a reduced Latin square. For j = 1, . . . ,n let Gj (resp. ttj) be 
the jth row (resp. column) of L. For any permutation a E Sn and any j = 1, . . . ,n, 
(a, avrjO"~ii^-,, avTj) G J^^^^ and any element of J^^j^ is of this form. 

Proof. Let 9 = (a,/3,7) be such that Q{L) G RLS{n). Since 7 is applied to all the 
symbols, the first row and first column of (q!,/3,1)(L) must be identical. Hence, after 
permuting the rows of L by a, /3 needs to permute the columns of (a, 1, 1){L) to yield a 
first row that is identical to one of the columns of (a, 1, 1)(-^)- Note that the first row 
of (a,l, 1)(L) is o"a-i(i) and the jth column of (a,l,l)(L) is anj (Lemma l(i)). Thus, 
after permuting the columns of (a, 1, 1)(-^) by /3 = aiXja'^^,-^. the first row of the resulting 
square is avr^ and, since avr^ is also one of the columns (it is was the jth column before 
the columns were permuted), it must be the first column of (a, /3, 1)(-^). Now, 7 needs to 
permute the symbols of anj to obtain the identity permutation (a reduced Latin square). 
By Lemma l(ii) aTij'')~^ = 1 and thus 7 = aiij. D 



When n is even, applying an isotopism preserves the parity of a Latin square 
Lemma 1), and thus the parity of all the Latin squares in an isotopy class is the same. 
When n is odd an isotopy class contains an equal number of even and odd Latin squares. 
This is not the case when only reduced Latin squares are concerned. In order to determine 
the number of even and odd reduced Latin squares in an isotopy class the following 
definition will be useful: 
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Definition 1. The parity type of a Latin square L of order n is defined to be {k,m), for 
^ k,m ^ n/2, if k (respectively m) of its rows (resp. columns) have one sign and the 
remaining n — k rows (resp. n — m columns) have the opposite sign. The parity class 
{k,m) consists of all Latin squares of parity type {k,m). 

Proposition 2. A parity class is the union of isotopy classes. 

Proof. Since an isotopy either changes the signs of all the rows (columns) or leaves them 
all unchanged, the parity type of two isotopic Latin squares is equal. D 

Thus, it is valid to refer to the parity type of an isotopy class, as parity type of any of 
its members. 

Lemma 2. Suppose n is odd. Let L E RLS{n) and let O = (a, aTTjCr^^i/-^-,, avTj) G J^n^, 
as in Proposition 1, then 

par(e(L)) = sgn(a«-i(i))sgn(7rj)par(L). 

Proof. By Proposition 3.1 in [10], par(6(L)) = (sgn(a))"(sgn(a;7rja~^i|.-^J)"par(L). Since 
n is odd, sgn(cr~ii^s) = sgn(ac,-i(i)) and sgn(a) appears an even number of times, the 
result follows. D 

The following theorem appears in a different form in [15] (Theorem 2.1). The proof 
here is different and was obtained independently. 

Theorem 1. Suppose n is odd and the parity type of L E RLS[n) is {k,m), then 



km+in-k)in^m) \^n^\ jf fc = m (mod 2) 
\^:,{L)nRELS{n)\ = { ^,^ J!^,^ J, :\^' _ ^ , \ _[ (2.1) 

■,2 



k{n-m)+m{n-k) \^i^(^^\ jf A;^m(mod2). 



Proof. Suppose k and m are either both even or both odd (first row in (2.1). Since n is 
odd, n — k and n — m are either both odd or both even respectively. Let I = {ii, . . . , i^} 
be the set of indices of the k rows of L having the same sign and let J = {ji, . . . ,jm} be 
the set of indices of the m columns of L having the same sign. Suppose L is even. Since 
k and m have the same parity, the rows indexed by / and the columns indexed by J must 
have the same sign and the rest of the rows and columns all have the opposite sign. Now, 
let = (a, ttTTjCr'ii^w avTj) e J^n,L (^^e Proposition 1) be such that 0(-^) is even. By 
Lemma 2, 0"q,-i(i) and ttj must have the same sign. So, either a~^(l) G / and j G J or 
a~"'^(l) G [?^] \ / and j G [n] \ J. Since the number of permutations a such that a""'^(l) G / 
(resp. a'^{l) ^ /) is {k/n)n\ = k{n — 1)! (resp. (n — k){n — 1)!), the number of isotopies 
G such that Q{L) is even is {km + {n — k){n — m)){n — l)\ and thus the number of distinct 
reduced even Latin squares that are isotopic to L is [km + {n — k){n — m)]{n — l)\/ |2l(L)|. 
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Hence, 

|^^,(L) n RELS{n)\ '^"^ + in - k){n - m)]{n - 1)! 



km 


+ {n 


2t(L) 
-k){n- 


m)]{n - 


-1)! 


n ■ n! 


km 


+ (n 


n ■ n\ 
-k){n- 


"^)] ^/ 


;(^)i 


|2l(L) 






n2 


^n 


5 



by (1.1). The other case is proved in an analogous manner and is left to the reader. D 

Remark 2. For \J^^{L) fl ROLS{n)\ the cases in (2.1) are reversed since /cm+ (n — k){n — 
777,) + k{n — m) + m,{n — k) = v? . 

Remark 3. Note that if /c = m = all the squares in J^ni^) ^^^ even. 

3 On even and odd Latin squares of odd order 

The following extension to all n of the Alon- Tarsi Latin square conjecture [1] appears in 

[15]: 

Conjecture 1. For all n, RELS{n) - ROLS{n) ^ 0. 

When n is even Conjecture 1 is equivalent to the Alon- Tarsi conjecture. When n is odd 
it is not known whether Conjecture 1 is equivalent to another extension of the Alon- Tarsi 
conjecture by Zappa [17]. 

Let Nn{k, m) be the number of reduced Latin squares of order n in the parity class (k, m). 
We have the following corollary of Theorem 1: 

Corollary 1. Suppose n is odd. The parity class [k^m] contains "^ 2 Nn{k^m) 

even reduced Latin squares if k = m (mod 2) and (."~"^)+"^(."~ ) Nn{k, m) even reduced 
Latin squares if k ^ m (mod 2) 

Proof. Since a parity class is the union of isotopy classes (Proposition 2), the result 
follows from Theorem 1 by summing over all the isotopy classes contained in a given 
parity class. D 

Corollary 1 yields an equivalent version of Conjecture 1 for odd n: 

Proposition 3. If n is odd Conjecture 1 is equivalent to 

Y^ (-l)'=+™(n-2fc)(n-2m)iV„(A;,m) 7^0. (3.1) 

k,m<n/2 
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Proof. By Corollary 1, if /c = m (mod 2) the difference between the number of even 
reduced Latin squares and the number of odd reduced Latin squares in the parity class 
{k,m) is 

km + (n — k)(n — m) kin — m) + min — k)\ ^^ ,^ . in — 2k)in — 2m) ,^ ,, , 

Nn[k,m) = Nn[k,m) 



n'^ n^ J 77,^ 

and when k ^ m (mod 2), this difference is — -^^^^ — ^"~ "^' Nn{k,m). Summing over all 
parity classes, 

RELSin) - ROLSin) = ^ V i-lf+'^in-2k)in-2m)NJk,m) (3.2) 

n ^-^ 

k,m<n/2 

Assuming Conjecture 1 statement (3.1) follows. D 

This calls for a study of the values Nni^k,m). When n = 5 only the parity classes 
(0,0) and (1, 1) are nonempty. When n = 6 the only empty classes are (0,2) and (2,0). 
When n = 7 all parity classes are nonempty, as is the case for n = 8, 9, 10 and 11 (Tables 
1-5). Thus, it seems reasonable to conjecture that for all n ^ 7 and all k,ni ^ n/2 
the parity class {k,m) is nonempty. Furthermore, an estimate for the values Nn{k,ni) 
is desired. It is conjectured in [3] that the distribution of the permutation that maps 
between two randomly chosen rows of a randomly chosen Latin square of order n, among 
all derangements, approaches uniformity as n — )■ oo (a derangement is a permutation 
without fixed points). Based on this, it seems reasonable to conjecture that for large 
enough n the distribution of the parities of the rows and columns of a random Latin square 
will be close to the random distribution. This should also hold when restricting to reduced 
Latin squares, namely, if the parities of the rows and columns were chosen randomly, the 
proportion of reduced Latin squares of parity class {k, m) among all reduced Latin squares 
would be rs(2) (^)/2^", where r = 2 (resp. s = 2) ii k < n/2 (resp. m < n/2) and 1 
otherwise, (since k,m ^ n/2). Thus, rs is always 4 if n is odd. 

Conjecture 2. If n is odd, then for all k,ni < n/2 

Nn{k,m) ^ ^^^\RLSin)l (3.3) 

and if k is even, then for all k,m ^ n/2 



"(DC) 

22 



N^{k,m) ~ -AhQ2nl\RLS{n)\ (3.4) 



where r = 1 (resp. s = 1) ii k = n/2 (resp. m = n/2) and 2 otherwise. 

Conjecture 2 is supported by the data in the following tables. Table 1 lists the real 
values of Nj{k,m), along with the expected values given by (3.3). The fourth column 
shows the ratios between the corresponding real and expected values. Tables 2-5 list the 
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{k,m) 


Njik,m) 


Expected Nj{k,m) 


N7[k,m} 


Expected N-j{k,m) 


(0,0) 


3270 


4136 


0.79 


(1,0) + (0,1) 


47040 


57904 


0.81 


(2,0) + (0,2) 


246960 


173720 


1.42 


(3,0) + (0,3) 


208740 


289536 


0.72 


(1,1) 


345450 


202676 


1.70 


(2,1) + (1,2) 


979020 


1216056 


0.81 


(3,1) + (1,3) 


2493120 


2026760 


1.23 


(2,2) 


1803690 


1824084 


0.99 


(3, 2) + (2, 3) 


5538960 


6080280 


0.91 


(3,3) 


5275830 


5066904 


1.04 



Table 1: Real and expected values of A^7(/c,m). 

counts of Nn{k,m), as calculated for 100,000 randomly selected reduced Latin squares 
(using the Jacobson-Matthews method [9]), for each of n = 8, 9, 10 and 11, along with the 
expected values calculated using (3.3) and (3.4) for this sample size (similar simulations 
were performed for n = 12, . . . , 20 and close values were observed as well for estimated 
and observed counts). The counts for the classes {k,m) and {m,k), when k y^ m, were 
added, since the corresponding classes have equal size. Note that as an estimated value 
does not mean that the class is expected to be empty, but that for the given sample size 
no member of the class is expected to be observed. 

Lemma 2.4 in [15] states that almost all Latin squares belong to parity classes {k, m) 
such that k,m^ n/63, and thus it can be viewed as a partial result for Conjecture 2. 

Proposition 4. Assuming Conjecture 2 holds then 

nodd \RLb[n)\ 



n— >oo 



Proof. Combining (3.2) and Conjecture 2 we have that for large odd n 
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{k,m) 



Ns{k,m) Expected Ns{k,m) 



Ns,{k,m) 



Expected Ng,{k,m) 



(0,0) 


11 


8 


1.38 


(0,1) + (1,0) 


94 


96 


0.98 


(0,2) + (2,0) 


344 


344 


1.00 


(0,4) + (4,0) 


421 


428 


0.98 


(0,3) + (3,0) 


705 


680 


1.04 


(1,1) 


394 


392 


1.01 


(1,2) + (2,1) 


2604 


2736 


0.95 


(1,4) + (4,1) 


3430 


3416 


1.00 


(1,3) + (3,1) 


5258 


5472 


0.96 


(2,2) 


5085 


4784 


1.06 


(2, 4) + (4, 2) 


12123 


11964 


1.01 


(4,4) 


7819 


7477 


1.05 


(2, 3) + (3, 2) 


19000 


19144 


0.99 


(3, 4) + (4, 3) 


23697 


23924 


0.99 


(3,3) 


17015 


19140 


0.99 



Table 2: Observed and expected values of Ng{k,m) for a sample of 100,000 randomly 
chosen reduced Latin squares. 



{k,m) 


N9{k,m) 


Expected Ng{k,m) 


NQ(k,m) 


Expected Ng{k,m) 


(0,0) 


2 





- 


(0,1) + (1,0) 


20 


24 


0.83 


(0,2) + (2,0) 


124 


112 


1.11 


(1,1) 


125 


124 


1.01 


(0,3) + (3,0) 


264 


256 


1.03 


(0,4) + (4,0) 


340 


384 


0.89 


(1,2) + (2,1) 


1031 


992 


1.04 


(1,3) + (3,1) 


2312 


2304 


1.00 


(1,4) + (4,1) 


3479 


3464 


1.00 


(2,2) 


1997 


1976 


1.01 


(2, 3) + (3, 2) 


9182 


9232 


0.99 


(2, 4) + (4, 2) 


13704 


13840 


0.99 


(3,3) 


10799 


10768 


1.00 


(3, 4) + (4, 3) 


32468 


32296 


1.01 


(4,4) 


24153 


24224 


1.00 



Table 3: Observed and expected values of Ng{k,m) for a sample of 100,000 randomly 
chosen reduced Latin squares. 
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{k, m) 



Nio{k,m) Expected A''io(A;, m) 



NiQ(k,m} 



Expected Nio{k,m) 



(0,0) 


2 


(0,1) + (1,0) 


8 


(0,2) + (2,0) 


32 


(1,1) 


40 


(0,3) + (3,0) 


86 


(0,5) + (5,0) 


107 


(0,4) + (4,0) 


165 


(1,2) + (2,1) 


335 


(1,3) + (3,1) 


881 


(1,5) + (5,1) 


939 


(2,2) 


769 


(1,4) + (4,1) 


1652 


(2, 3) + (3, 2) 


4060 


(2, 5) + (5, 2) 


4181 


(2, 4) + (4, 2) 


7326 


(3,3) 


5400 


(3, 5) + (5, 3) 


11461 


(5,5) 


6083 


(3, 4) + (4, 3) 


19208 


(4, 5) + (5, 4) 


20286 


(4,4) 


16979 




8 

32 

40 

88 

96 

160 

344 

912 

960 

772 

1600 

4120 

4324 

7208 

5492 

11536 

6056 

19224 

20188 

16824 



1.00 
1.00 
1.00 
0.98 
1.11 
1.03 
0.97 
0.97 
0.98 
1.00 
1.03 
0.99 
0.97 
1.02 
0.98 
0.99 
1.00 
1.00 
1.00 
1.01 



Table 4: Observed and expected values of NiQ{k,'m) for a sample of 100,000 randomly 
chosen reduced Latin squares. 



Now, 



k,m=0 



L'V ^J / \ / 

^(-l)^+-(n-2fc)(n-2m)r r 

. m—n \ / \ 

n/2j , . 

Y.{-l)\n-2k)r\ 

I.— n ^ ^ 



(3.6) 



fc=0 

\n/2\ 



n 



E(-i) 



A;=0 ^ ^ fc=0 

It can be shown by induction on n that for odd n ^ 3 

\n/2\ 



L"/2J / X 

2 E (-!)'<:) 



E(-i)'*(l) = (-i)^"^""' 



fc=0 



n-1 
2 
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{k, m) 



Nii{k,m) Expected A''ii(A;, m) 



Nii(k,mf 



Expected Nii{k,m) 



(0,0) 





(0,1) + (1,0) 


3 


(0,2) + (2,0) 


13 


(0,3) + (3,0) 


40 


(1,1) 


11 


(0,4) + (4,0) 


59 


(0,5) + (5,0) 


90 


(1,2) + (2,1) 


120 


(1,3) + (3,1) 


331 


(2,2) 


293 


(1,4) + (4,1) 


691 


(1,5) + (5,1) 


985 


(2, 3) + (3, 2) 


1699 


(2, 4) + (4, 2) 


3312 


(2, 5) + (5, 2) 


4881 


(3,3) 


2555 


(3, 4) + (4, 3) 


10209 


(3, 5) + (5, 3) 


14446 


(4,4) 


10472 


(4, 5) + (5, 4) 


29244 


(5,5) 


20546 





8 

32 

12 

64 

88 

112 

344 

288 

696 

968 

1728 

3464 

4848 

2596 

10385 

14536 

10384 

29080 

20356 



1.63 
1.25 
0.92 
0.92 

1.02 
1.07 
0.96 
1.02 
0.99 
1.02 
0.98 
0.96 
1.01 
0.98 
0.98 
0.99 
1.01 
1.01 
1.01 



Table 5: Observed and expected values of Nii{k,'m) for a sample of 100,000 randomly 
chosen reduced Latin squares. 



and 



\n/2\ 

El 

fc=0 



n 



n-iin-2 
-1^ ' I n-l 



Thus, the sum in (3.6) is and the result follows. 



D 



The conjecture presented in Proposition 4 is supported by Theorem 3.2 in [15] which 
states that \RELS{n)\ and \ROLS{n)\ share a large divisor which grows exponentially as 



n — )■ oo. 



4 On the order of 2t(L) 

Suppose L and L' are isotopic Latin squares. Note that if L' = Q{L), then for any F G J^n, 
T G 2l(L) if and only if 616-1 G 2l(L'). Hence |2l(L)| = |2l(L')|. Since every Latin square 
is isotopic to some reduced Latin square, it is enough to study the orders of the autotopy 
groups of reduced Latin squares. 
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The parity type of a square provides some information about its autotopy group. Brown- 
ing, Stones and Wanless ([2], Lemma 4.1) showed that for a prime order p there is exactly 
one isotopism class that contains Latin squares L for which p^ divides |2l(L)|, namely, the 
isotopism class containing the Cayley table of Zp. The following corollary of Theorem 1 
extends that result: 

Corollary 2. Let L be a Latin square of prime order p > 2 and of parity type {k,m). 

(i) If k ^0 orm y^O then p^ /|2t(L)|. 

(ii) If k and m are both nonzero then p /|2l(L)|. 

(iii) // the parity of L is (0,0) but L is not isotopic to the Cayley table of Zp, then 

p' j(mL)\ 

Proof. Statements (i) and (iii) follow from the above mentioned result in [2] and the fact 
that the parity type of a Cayley table of odd order is (0, 0). Without loss of generality we 
may assume that L is reduced (otherwise we take a reduced Latin square in the isotopy 
class of L and use the fact that isotopy preserves the parity type). Setting n = p we 
look at one of the expressions for ^^^{L) fl RELS(p)\ in (2.1). If just one of k and m is 
nonzero, then p divides the numerator but p^ does not. Thus, p must divide | J^p(L)| and 
by (1.1) the highest power oi p that may divide |2l(L)| is 1. This proves (ii). D 

For the proof of the next theorem the following definition will be useful: 

Definition 2. Let I, J G [n] with |/|, \J\ ^ n/2. We say that a Latin square L of order 
n has parity set (/, J) if the rows (resp. columns) indexed by / (resp. J) are all of the 
same parity and the rows indexed by [n] \ I (resp. [n] \ J) are all of the opposite parity. 

Remark 4. Permuting the columns of a Latin square preserves the first component of its 
parity set and permuting the rows of a Latin square preserves the second component of 
its parity set. Permuting the entries of a Latin square preserves its parity set. 

Theorem 5. Let L be a Latin square of order n with parity type {k,m), then 

// no column of L is the product of two other columns, then 

mL)\ ^ --. (4.2) 

[Similarly, if no row of L is the product of two other rows, then |2l(L)| ^ n!/(^).) 
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Proof. We can assume that L is a reduced Latin square. By Remark 4, choosing all 
different a G S'„, in (a, a7^jO"~il^^, avTj) G '^n,L (^^^ Proposition 1) will yield all (^) 
distinct / C [n] of order k as the first component of the parity set of the reduced Latin 
squares that are isotopic to L. Thus |^^(L)| ^ Q. By (1.1) |2l(L)| ^ ^. Since 

the same argument applies by looking at the second component of the parity set, the 
inequality (4.1) follows. Let S = {ai, . . . ,a(n\} be a set of permutations that, when 

taken as the first component of corresponding (^) isotopies on L, will yield (^) distinct 
reduced Latin squares with distinct first parity set component. Now, let a G S and 



suppose that 6i = (a, a7rjcr^_i^w avTj) G J^n,L ^^^ ©2 = {oi,a7iia^ 
satisfy ei(L) = e2(L). Let $ = e^^Gi Thus $ G 2l(L). That is. 



-HI) 



$ 



.1,0-0 



Ki)^i 



J o 



-1(1) 



TT, TTi 



G 2t(L). 



«7r;) G y;^L 



(4.3) 



For simplicity, denote a = crQ,-i(i). Let /3 = an^^Tija and 7 = tt^^t^j = <J~^f3a. We have 
that $ = (1, /3, cr~^/3cr) G 2l(L). While applying $ to L, after permuting the columns by 
P, the first row is /3 (by Lemma l(i), since it was originally the identity permutation) and 
we have to apply again /3 on the symbols in order to transform the first row back into 
the identity permutation (Lemma l(ii)). Thus 7 = /3 and we have $ = (1, (3, (3). Suppose 
/3(1) = r. After permuting the columns by /3, the original first column (the identity 
permutation) moves to the rth position. Then when applying /3 on the symbols The rth 
column becomes /3~^ (Lemma l(ii)), and since we assumed that (l,/9,/3) G 2l(L) we must 



have that (3 



TTr 



But (3- 



r 



71 ■ Til. Thus TTi 



TTjTTj. 



If we assume that no 



column of L is the product of two other columns, we must have that 6i(-^) 7^ 02 (-^). We 
conclude that for any a G S the different n isotopies obtained by taking the n columns 
of L yield n distinct squares that are isotopic to L. Since for different elements of S the 
squares obtained have parity sets with distinct first component, \J^^{L)\ ^ n • (^). By 

(i.i)|a(L)K^ = ^ 

Example 1. The parity type of the following Latin square of order 6 is (3, 3) and thus, 
by (4.1), |2l(L)| ^ (6 ■ 6!)/(3) = 216. Indeed, the autotopy group of this square has size 
216, so it is possible to reach the bound in (4.1). 



1 


2 


3 


4 


5 


6 


2 


1 


4 


3 


6 


5 
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5 
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1 



5 Computing 2t(L) by cycle structures 

Every permutation a E Sn can be decomposed into a unique (up to order) product of 
disjoint cycles. Cycle structures of permutations were considered in the context of Latin 
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squares in different aspects. Cavenagh, Greenhill and Wanless [3] considered the cycle 
structure of the permutation that transforms one row of a Latin square to another row. 
Other works [7, 14] considered the cycle structure of the permutations a, (3 and 7 in an 
isotopism O = (a,/3,7), in order to derive information on Latin squares for which G is 
an autotopism. Here the cycle structure of the rows of a Latin square L are considered 
in order to derive information on 2l(L). 

Lemma 3. Let L be a reduced Latin square of order n with rows {(Ti}"^^. Let = 
(a,/3,7) G 2l(L). Then for each i = 1, . . . ,n, 0"~-i(i-)CrQ-i(j) has the same cycle structure 
as ai. 

Proof. By Proposition 1, has the form = (a, Q^^jcr^-i(i); ci'^j) ^ 2t(L) for some column 
TTj of L. We follow the zth row of L as is applied. Originally it is ai. After permuting 
the rows the zth row is ^^-'^-(i) (following the convention that i in the jth place of a 
permutation a signifies that a{i) = j). After permuting the columns by aTija~\,^. the 
ith row becomes aTija~\,^.aa-^{i) (Lemma l(i)) and after permuting the symbols, the 

resulting ith. row is a7^jC^~il^^(Tc^-l(^)(«v^j)~^ (Lemma l(ii)). Since is an autotopism of 
L we have 

ai = a7rja~l^^^^aa-i(i)ianjy^ (5.1) 

Thus ai and cra-^n)0'a-^(i) are conjugates and hence have the same cycle structure. D 

Theorem 6. The following is an algorithm for finding 2t(L) for a given reduced Latin 
square L of order n. 

(1) Compute the cycle structures of the rows of L, viewed as permutations in Sn- Let 
Ci, . . . ,Cs be the distinct cycle structures of the rows, sorted in some well-defined 
way (see [7]) and let A = (Ai, A2, . . . , A^) be a partition of n where each Aj is the 
number of rows with cycle structure Ci. 

(2) For each row au, compute the cycle structures of the set of permutations {a^^ai}^^]^. 
Let Ci, . . . , Cf be the cycle structures of these permutations, sorted as in (1), and 
let \^ be the partition of n corresponding to these cycle structures. 

(3) // the ordered sets {Ci, . . . , Cg} and {Cf , . . . , C^} coincide and X = X^ construct the 
set Ik of permutations a & Sn satisfying: 

(*) For each i, aa{i) and a^^ai have the same cycle structure. 

(4) For each a G Ufc/^ and each column iij, j = 1, . . . ,n, let Qaj = («? ^^j'^Q-inv ^^j)- 
Check whether Qa,j{L) = L. If equality holds then Qaj G 2l(L). 

Proof. The condition (*) follows directly from Lemma 3. Since a is a bijection the con- 
dition in step (3) must hold. D 

Remark 7. Condition (*) implies that a{k) = 1 since the identity permutation has its 
own unique cycle structure. 
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Example 2. Consider the following reduced Latin square L of order 7: 
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The cycle representations of the rows of L, grouped by cycle structure, are: 

Row Cycle representation 



(5.2) 



1 (1)(2)(3)(4)(5)(6)(7) 

2 (1,2)(3,4)(5,7,6) 

4 (1,4)(3,7)(2,5,6) 

5 (1,5)(3,6)(2,4,7) 
7 (4,6)(1,3,5,2,7) 

3 (1,6,7,5,4,2,3) 

6 (1,7,4,5,3,2,6) 

One set of possible a consists of permutations that satisfy oi{l) = 1, a(7) = 7 and permute 
the sets {2,4,5} and {3,6}. This set of permutations contains 3! ■ 2! = 12 permutations. 
(This corresponds to A; = 1 in Step (2) of the algorithm in Theorem 6.) 
If we take the product of each row by a^^^ (the inverse of the 4th row, /c = 4 in Step (2) 
of the algorithm in Theorem 6) we obtain the same set of cycle structures: 

i Cycle rep. of cr^'^ai 



(5.3) 



4 (1)(2)(3)(4)(5)(6)(7) 


1 (1,4)(3,7)(2,6,5) 
3 (1,5)(2,7)(3,4,6) 
5 (1,2)(6,7)(3,5,4) 


7 (2,3)(1,7,4,5,6) 


2 (1,6,2,4,7,5,3) 
6 (1,3,6,4,2,5,7) 



So, another set of possible a's consists of permutations that satisfy a(4) = 1, a(7) = 7 and 
map the set {1,3,5} onto the set {2,4,5} and the set {2,6} onto the set {3,6}. Again, 
there are 3! ■ 2! = 12 such permutations. Taking the product of each row by a^^ also 
produces the same set of cycle structures as the original one, while taking the products 
of the rows with each of a^^, for A; = 2, 3, 6, 7, produces sets of cycle structures that do 
not coincide with the original one. Thus, |Ufc/fc| = 3 ■ 12 = 36 and there are 36 ■ 7 = 252 
isotopisms to check in step 4, of which only three are autotopisms. Hence |2l(L)| = 3. 



Remark 8. The algorithm, as described in Theorem 6, was tested on 100,000 randomly 
selected Latin squares (Jacobson-Matthews method [9]) of each of the orders 8,9,10, and 
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11. It performed faster than the algorithm by McKay, Meynert and Myrvold [11], based 
on "nauty" , but slower than an improved version of "nauty" using vertex invariants (such 
as the "train" [16]). 

The algorithm of Theorem 6 can be massively sped up by constructing the permuta- 
tions a in Step (4), in parts, corresponding to the cycle structures Cf in Step (3). After 
constructing each part we can perform Step (4) on the rows of L that are permuted by 
the part of a already constructed. This may rule out most of the candidates, or produce 
candidates a while saving the time needed to consider all the parts corresponding to the 
cycle structures C^. This idea is illustrated in the following example. 

Example 3. Consider the Latin square in Example 2. We saw that if we take the product 
of each row with crj^ the same set of cycle structures as the original one is obtained (see 
(5.2) and (5.3)). This implies that one set of possible a may satisfy a(7) = 7. By (5.1) it 
follows that 

a-j = aiija'^^a-jTxJ^a"^ (5.4) 

for one or more columns tvj. Let r = a'^^a-j = (2, 3)(1, 7, 4, 5, 6) (cycle decomposition). 
Since tti = 1 (the identity permutation) we have ttittt]"^ = r and by (5.4) we have 

aj = (4, 6)(1, 3, 5, 2, 7) = (a(2)a(3))(a(l), a(7), a(4), a(5), a(6)). (5.5) 

Since we know that a(4) = 1 and a{7) = 7, by looking at the second cycle in (5.5), it 
follows that a{5) = 3 and a(6) = 5, both in contradiction to the constrains in Example 2 
requiring that a must map 5 to one of {2, 4, 5} and 6 to one of {3,6}. Thus all possible 
a for /c = 4 and j = 1 in Steps (2)-(4) of Theorem 6 are ruled out. 

If we take j = 2 in Step (4) of Theorem 6 we obtain 7r2r-7r^^ = (1, 6)(2, 5, 3, 7, 4) and 
thus, by (5.1), it follows that 

ar = (4, 6)(1, 3, 5, 2, 7) = («(1), «(6))(«(2), «(5), «(3), a(7), «(4)). (5.6) 

Using the fact that a(4) = 1 and a(7) = 7 it follows that a = (1,4)(2,3) is the only 
permutation that coincides with the constrains in Example 2. It remains to check whether 
the isotopy Qa,2iL) = {a, a7i20'~}.i,^-., a7i2) satisfies Qa,2{L) = L. 

The reader can verify that taking the row n^ produces another candidate a = 
(1, 5,4)(2, 6, 3), while taking the rows 713,71^,71^ and 71-1 yield contradictions. Thus, for 
the case fc = 4 in Step (2), the list of possible autotopisms to check in Step (4) was 
narrowed down to only 2 in n = 7 iterations instead of ^t- 11^=1 AJ = 7 ■ 2! ■ 3! = 84 in 
Y[i=i -^J = 2! ■ 3! = 12 iterations according to the original algorithm. 

Corollary 3. Let L be a Latin square of order n. Let Ci, . . . ,Cs be the distinct cycle 
structures of the rows and let A = (Ai, A2, . . . , A^) be a partition of n where each Aj is the 
number of rows with cycle structure Ci . Then 

s 

mL)\^n'l[Xi\. (5.7) 
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// no column of L is the product of two other columns then 



|a(L)|^nJ]A,! 



(5. 



j=i 



Proof. Without loss of generality we may assume that L is reduced. For each k satisfying 
the condition of step (3) in Theorem 6, there are ni=i -^i' permutations a satisfying (*). 
Since there are n rows to check in step (2) and for each a there are n isotopies Q^j 
to check in step (4) of Theorem 6, Inequality 5.7 follows. It was shown in the proof of 
Theorem 5 that if two distinct isotpisms, with the same first component a, map L to the 
same reduced Latin square, then L must have a column that is the product of two other 
columns. If no such column exists, then in step (4) of Theorem 6 at most one value of j 
satisfies 6qj(L) = L. Thus Inequality 5.8 follows. D 

Remark 9. The worst case in the algorithm in Theorem 6 is when n ~ 1 rows have the 
same cycle structure (the first row, the identity permutation, has its own cycle structure). 
This condition is achieved, for example, by atomic Latin squares [16] and by Cayley tables 
of elementary abelian groups. In this case the bound in (5.7) is n • n!, which provides no 
information, as {J^n l\ = n ■ n\. 

Example 4. The rows of the Latin square 
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have cycle structures 



Row Cycle representation 

^~ (1)(2)(3)(4)(5)(6) 

2 (1,2)(3,4)(5,6) 

3 (1,5,3)(2,6,4) 

5 (1,3,5)(2,4,6) 

4 (1,6,3,2,5,4) 

6 (1,4,5,2,3,6) 

The products of the rows with each o"^^, /c = 1, . . . , 6, produce a set of permutations with 
the same cycle structures as above. Thus, the set U^/fc has size 2! ■ 2! ■ 6 = 24, and the 
bound in (5.7) is 24 ■ 6 = 144. The actual size of 2l(L) in this case is 72. This is the 
closest example found, among Latin squares of order 6 and 7, to the bound in (5.7). 
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